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ABSTRACT 

We present an observational and dynamical study of newly discovered main-belt comet 313P/Gibbs. 
We find that the object is clearly active both in observations obtained in 2014 and in precovery 
observations obtained in 2003 by the Sloan Digital Sky Survey, strongly suggesting that its activity 
is sublimation-driven. This conclusion is supported by a photometric analysis showing an increase in 
the total brightness of the comet over the 2014 observing period, and dust modeling results showing 
that the dust emission persists over at least three months during both active periods, where we find 
start dates for emission no later than 2003 July 24T10 for the 2003 active period and 2014 July 28T10 
for the 2014 active period. Erom serendipitous observations by the Subaru Telescope in 2004 when 
the object was apparently inactive, we estimate that the nucleus has an absolute R-band magnitude 
of Hr = 17.1 ± 0.3, corresponding to an effective nucleus radius of Tg ^ 1.00 ± 0.15 km. The object’s 
faintness at that time means we cannot rule out the presence of activity, and so this computed radius 
should be considered an upper limit. We find that 313P’s orbit is intrinsically chaotic, having a 
Lyapunov time of Ti = 12 000 yr and being located near two 3-body mean-motion resonances with 
Jupiter and Saturn, 11J-1S-5A and 10J+12S-7A, yet appears stable over >50 Myr in an apparent 
example of stable chaos. We furthermore find that 313P is the second main-belt comet, after P/2012 
Tl (PANSTARRS), to belong to the ^155 Myr old Lixiaohua asteroid family. 

Subject headings: comets: general — comets: individual (313P/Gibbs) — minor planets, asteroids: 
general 


1. INTRODUCTION 

Comet 313P/Gibbs was discovered as P/2014 S4 
(Gibbs) on 2014 September 24 by A. R. Gibbs using the 
Catalina Sky Survey’s 0.68-m Schmidt telescope. Eollow- 
up images from the Mt. Lemmon 1.5 m telescope revealed 
a 20" tail and a slightly elliptical coma with a full-width 
at half-maximum (EWHM) approxima tely twice that of 
nearby field stars (Gibbs & Sato [2014 ). Analysis of ob¬ 
servations by the Hubble Space Telescope (HST), Keck 
I, and the Danish 1.5 m telescope at La Silla indicate 
continuous mass loss between 2014 October 2 and 2014 
November 6, inconsistent with an impact event, though 
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no spectroscopic e vidence of gas emission was detected 
(jJewitt et al.|2015|). As of 2015 January 15, 313P has an 
osculating semimajor axis, eccentricity, and inclination 
of a = 3.156 AU, e = 0.242, and i = 10.97°, respectively, 
an orbital period of 5.61 yr, and a Tisserand parameter 
with respect to Jupiter of Tj = 3.132. These orbital el¬ 
ements place the object unequivocally within the main 
asteroid belt, maki ng it one of the growing number of 
active asteroids (cf. |Jewi in the main belt to be 


discovered in recent years. 

Active asteroids inc lude main-belt comets (MBGs; cf. 


Hsieh & Jewitt 2006), which exhibit comet ary a ctivity 
due to ice su blirnation, and disrupted asteroids (cf. [Hsieh 


et al. 2012a), where dust emission is caused by ettects 


such as impacts or rotational disruptions. Distinguishing 
between different sources of activity is crucial for devel¬ 
oping our understanding of the global properties (e.g., 
abundances, spatial distributions, and physical charac¬ 
teristics) of objects exhibiting each type of activity (e.g., 
Hsieh etal.jl^ . 


Unambiguously determining the source of comet-like 
activity, especially for a recently discovered object, is 
unfortunately not straightforward. Des pite detections of 


water ice frost on asteroid ( 24) Themis (Rivkin & Emery 


2010 Campins e t al.|2010) and water vapor outgassmg 
from (1) Ceres ( jKiippers et al.lpMil, attempts to di- 


troscopy (e.g., Jewitt et al. |2009| Licandro et al. 

2011 

Hsieh et al. |2012a|b|cl 

20131 |de Val-Borro et al. 

2012 

O’Rourke et al. 2013|) 

lave all been unsuccessfu 

. We 
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photometric monitoring (cf. |Hsieh et al. 2012a), to de¬ 
termine whether dust emission w as longdived (imply ing 
sublimation-driven emission; e.g., Hsieh et al.||2Q12b| or 
■pulsive (imply ing impact-driven ejection; e.g., ( 


teven- 


son et al. 2012|). As such, unambiguous confirmation 
ot sublimation is difficult to achieve, a situation further 
complicated in ca ses where multiple m echanisms could 


be at work (e.g., Jewitt et al. 
tivity near periheli on, such as 


2014a|). Recurrent ac- 
lat observed for MBCs 


2004 2010 


Jewitt et al. 


133P/ Elst-Pizarro (Hsie h et al. _ 

2014bj) and 238P/Read (|Hsieh et al.||2011bpj is currently 
considered the most reliable indicator of sublimation- 
driven activity. This behavior is naturally explained by 
thermally-modulated sublimation and extremely difficult 
to explain as a consequence of any physical disruption 
( Hsieh et al.||2012a Jewitt|[2M2 ). 


2. OBSERVATIONS 

Following 313P’s discovery, we identified pre-discovery 
images from August 2014 obtained by the 1.8 m Pan- 
STARRSl (PS 1) survey telescope using PSl ipi- and 
rcpi-band filters (Tonry et al. 2012). We also obtained 
follow-up imaging in October and November 2014 us¬ 
ing the Large Format Camera (Simcoe et al. 2000) on 
Palomar Observato ry’s 5.1 m Hale telescope, the Large 
Monolithic Imager (Bida et al.|2014) on Lowell Observa¬ 
tory’s 4.3-m Discovery Channel T'elescope, a 4096 x 4096 
Fairchild CCD on the 2.0 m Faulkes Telescope North on 
Haleakala, a 2048 x 2048 Textronix CCD on the Univer¬ 
sity of Hawaii (UH) 2.2 m telescope on Manna Kea, and 
the European Southern Obs ervatory (ESQ) Faint Object 
Spectrograph and Camera (Buzzoni et al. 1984|) on the 
3.58 m New Technology Telescope (Program 194.C-0207) 
operated by ESO at La Silla. SDSS r'-band or Kron- 
Cousins R-band filters, as specified in Tablewere used 
for all follow-up observations. We used Imag e Reduction 
and Analysis Facility software (|Tody 1986|) to perform 
standard bias subtraction and flat-held reduction (using 
dithered twilight sky images) for all data, except those 
from PSl, which wer e reduced usin g the system’s Image 
Processing Pipeline ( Magnier||2006 ). 

Using the Solar System Object Image Search tool 
(SSOIS; |Gwyn et al. 2012), provided by the Canadian 
Astronomical Data Centre, we identified precovery ob¬ 
servations obtained in 2003 b y the 2.5 m Sloan Digital 
Sky Survey (SPSS) telescope (York et^Mg^Ql iFnkngita 


et al.||1996| |Gunn et al.||1998| 2006| lAihara et al.pO H )" 

and in 2004 by Suprime-Carn (jlVliyazaki et al. 2002 1) on 
the 8.2 m Subaru Telescope. The SDSS observations, ob¬ 
tained in g'-^ r'-, i'-^ and 2 ;'-band Alters, show a clearly 
comet ary object with a dust tail extending as much as 
^1 arcmin in the antisolar direction. The object is ex¬ 
tremely faint in the Subaru images, obtained using a 
wide-band “W-J-VR” Alter encompassing both the John¬ 
son V and Cousins R passbands, exhibiting no indica¬ 
tions of activity. 

Absolute calibration of all data was accomplished using 
magnitudes of held stars with ap proximat ely solar colors 
from SDSS (lAihara et ffil|2Qll|) or PSl dSch lafly et al. 
i Tonry et al. 2U12[ Magnier et al. 2U13| . Conver- 


2012 


Sion of all photometry to R- band was accomplish ed using 
transformations derived by |Tonry et al. (2012) and by 
R. Lupton (http://www. sdss.org/). Comet photome¬ 
try was performed using circular apertures with varying 


radii depending on the nightly seeing, where background 
statistics were measured in nearby regions of blank sky, 
avoiding contamination from the comet. 

3. RESULTS 

3.1. Photometric Analysis 

Apparent R-band magnitudes as measured in 5" cir¬ 
cular apertures, and equivalent absolute magnitudes (as¬ 
suming a standard iL, G phase function where we assume 
G = 0.15), as well as absolute magnitudes correspond¬ 
ing to larger rectangular apertures encompassing the en¬ 
tire comet, for all observations are listed in Table 
Both the central and total brightnesses of the comet 
decline over the 2003 observing period (which spans a 
true anomaly range of 30.5° < n < 46.7°), declining 
by ^40% in 59 days. Meanwhile, the comet brightens 
over the 2014 observing period (—6.5° < v < 24.0°), in¬ 
creasing in total brightness by ^150% in 104 days (Fig¬ 
ure IT). The object’s increasing brightness in 2014 is 
suggestive of ongoing, sublimation-driven dust produc- 
tion, but is not u nambiguous proof of sublimation (cf. 
Hsieh et al.||2012a|). Multi-filter SDSS observations from 
2003 show that the comet has approximately solar col¬ 
ors {g' -r' = 0.54 ± 0.05, r' - i' = 0.10 ± 0.05, and 
i' — z' = 0.05 ± 0.05; equivalent to B — V = 0.71 ± 0.05, 
V -R = 0.41 ± 0.05, and R - / = 0.33 ± 0.05) with no 
appreciable color differences between the coma and the 
tail, suggesting that the measured colors in both regions 
are dominated by dust. Strong coma during both ac¬ 
tive periods prevent any meaningful determination of a 
rotational period or amplitude. 

We measure an apparent R-band magnitude of tur = 
23.3 ± 0.2 for 313P in 2004, when it was apparently in¬ 
active, while it was at u = 108.7° and at a heliocentric 
distance of R = 3.21 AU. This corresponds to an ab¬ 
solute magnitude of Hr = 17.1 ± 0.2, from which we 
compute an effective nucleus radius of ^ 1.0 ± 0.1 km 


(assumin g a R-band albedo of pr = 0.05; cf. Hsieh et 
al.||2009l), making it a mid-sized MBG (cf. |Hsieh et al 
hile no coma or tail was visible in this data 


the object’s faintness at the time means that we cannot 
rule out the presence of activity, particularly since ac¬ 
tivity was detected for MBGs 133P and P/2010 R2 (La 
Sagra) as late as = 109° and n = 117°, respectively. 
when both were at R ^ 3.25 AU (Kaluna et al. [2011 ' 


Hsieh||2014|). In fact, HST observations in dicate 313P s 
nucle us may be as small as Vq ^ 0.5 km (| Jewitt et al. 
2015). As such, the nucleus size calculatea here should 
be regarded as an upper limit. 

Using this nucleus size, we can estiniate to tal dust 
masses at different times. Following |Hsieh| ( [2014| ), we And 
a peak total dust mass of ^ (2.5 ±1.0) x 10^ kg and 
dust-to-nucleus mass ratio of M^/Mn ~ (4.5±1.5) x 10“^ 
in 2003 (at u = 30.5°) and Md ^ {3 ± 1) x 10^ kg and 
MdiMN - (5.5 ± 2.0) X 10-^ in 2014 (at z/ = 24°), or 
a decline in peak total dust mass of ^ 1 order of mag¬ 
nitude between 2003 and 2014. Additional photometric 
measurements through the end of the 2014-2015 observ¬ 
ing window will permit a more direct comparison of 2003 
and 2014-2015 activity levels over the same orbital arcs, 
while more accurate measurements of 313P’s nucleus size 
will improve the accuracy of the dust mass estimates re¬ 
ported here. However, these current results suggest that 
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Fig. 1.— Orbital position plot (left) and composite images (right) constructed from observations detailed in Tablewith the 
Sun (black dot) at the center of the orbital position plot, and the orbits of Mercury, Venus, Earth, Mars, 313P, and Jupiter 
shown as black lines. Perihelion (P) and aphelion (A) are marked with crosses. Colored points (left panel) and images (right 
panel) correspond to observations from (a) 2003 September 30, (b) 2003 October 23, (c) 2003 November 28, (d) 2014 August 6, 
(e) 2014 September 3, (f) 2014 October 1, (g) 2014 October 30, and (h) 2014 November 18. An open square symbol (left panel) 
marks the position of observations from 2004 September 16 (image not shown). The object is at the center of each image panel, 
and North (N), East (E), the antisolar vector (—O), and the negative heliocentric velocity vector (—r’) are marked with arrows. 


TABLE 1 
Observations 


UT Date 

Tel.^ 

Nb 

t° 

Filter 

R^ 

A® 



(R, A, a)^ 


rriR^toti^: 1,0)J 

Active?^ 

2003 June 20 

Perihelion... 



2.367 

2.604 

23.0 

0.0 

— 

— 

— 

— 

2003 September 30 

SDSS 

4 

216 

g'r'i'z' 

2.434 

1.598 

16.1 

30.5 

18.4T0.1 

14.6T0.1 

13.6T0.1 

yes 

2003 October 23 

SDSS 

4 

216 

g'r'i'z' 

2.465 

1.512 

8.5 

37.0 

18.2T0.1 

14.8T0.1 

13.9T0.1 

yes 

2003 October 24 

SDSS 

4 

216 

g'r'i'z' 

2.467 

1.511 

8.3 

37.2 

18.3T0.1 

14.9T0.1 

14.1T0.1 

yes 

2003 November 28 

SDSS 

4 

216 

g'r'i'z' 

2.525 

1.636 

12.1 

46.7 

18.7T0.1 

14.9T0.1 

14.2T0.1 

yes 

2004 September 16 

Subaru 

2 

120 

W-J-VR 

3.208 

3.742 

14.1 

108.7 

23.3T0.2 

17.1T0.2 

— 

no 

2009 January 18 

Perihelion... 



2.385 

2.977 

17.0 

0.0 

— 

— 

— 

— 

2014 August 29 

Perihelion... 



2.392 

1.555 

16.9 

0.0 

— 

— 

— 

— 

2014 August 06 

PSl 

2 

90 

WY>1 

2.395 

1.750 

22.1 

353.5 

20.8T0.1 

16.6T0.1 

16.6T0.1 

yes 

2014 August 18 

PSl 

3 

135 

Zpi 

2.392 

1.637 

19.6 

357.0 

20.6T0.1 

16.6T0.1 

16.6T0.1 

yes 

2014 September 03 

PSl 

2 

90 

WY>1 

2.392 

1.518 

15.2 

1.8 

20.3T0.1 

16.7T0.1 

16.7T0.1 

yes 

2014 October 01 

P200 

38 

1140 

r' 

2.399 

1.429 

7.7 

10.0 

19.3T0.1 

16.1T0.1 

15.8T0.1 

yes 

2014 October 13 

ETN 

2 

120 

r' 

2.404 

1.444 

8.4 

13.2 

19.4T0.1 

16.1T0.1 

16.0T0.1 

yes 

2014 October 17 

DCT 

29 

2610 

r' 

2.407 

1.461 

9.6 

14.7 

19.3T0.1 

15.9T0.1 

15.6T0.1 

yes 

2014 October 22 

DCT 

2 

240 

R 

2.410 

1.484 

11.0 

16.2 

19.6T0.1 

16.1T0.1 

15.8T0.1 

yes 

2014 October 30 

UH2.2 

5 

1500 

R 

2.416 

1.530 

13.4 

18.5 

19.6T0.1 

16.0T0.1 

15.8T0.1 

yes 

2014 November 17 

NTT 

4 

480 

r' 

2.432 

1.683 

18.4 

23.7 

20.2T0.1 

16.2T0.1 

16.0T0.1 

yes 

2014 November 18 

NTT 

3 

900 

r' 

2.433 

1.694 

18.6 

24.0 

20.6T0.1 

16.6T0.1 

15.7T0.1 

yes 

2014 December 16 

UH2.2 



R 

2.464 

2.023 

22.7 

31.9 




yes 

2020 April 15 

Perihelion... 



2.419 

3.344 

7.8 

0.0 

— 

— 

— 

— 


^ Telescope (SPSS: Sloan Digital Sky Survey; PSl: Pan-STARRSl; P200: Palomar Hale Telescope; FTN: Faulkes Telescope North; DCT: 
Discovery Channel Telescope; UH2.2: UH 2.2 m telescope; NTT: New Technology Telescope). 

^ Number of exposures. 

Total integration time, in s. 

^ Heliocentric distance, in AU. 

® Geocentric distance, in AU. 

^ Solar phase angle (Sun-object-Earth), in degrees, 
s True anomaly, in degrees. 

^ Mean apparent R-band magnitude, assuming solar colors. 

^ Absolute R-band magnitude (at R = A = 1 AU and 0 = 0°) as measured with a 5'.'0 aperture, assuming lAU H^G phase-darkening 
where G = 0.15. 

j Absolute R-band magnitude as measured with rectangular aperture enclosing entire visible extent of comet. 

^ Is visible activity present? 
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Fig. 2.— (a) Plot of absolute total R-band magnitudes measured for 313P (including the total observed flux from the nucleus, 
coma, and dust tail) from 2003 through 2014. (b) Predicted apparent R-band magnitude of the nucleus of 313P through the 
end of the 2016-2017 observing window. Shaded regions mark time periods when the object is not observable from Earth (i.e., 
solar elongation <50°). 
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there is a significant decline in activity strength after 
just two orbit passages, perhaps due to a combination of 
mantling of the active area re sponsible for the observed 
activity (cf. Hsieh et 'ar]|2Q15| ) and depletion of exposed 
volatile material. Alternatively, other effects (e.g., out¬ 
bursts or pole orientation evolution) could be responsi¬ 
ble, and as such, monitoring of future active episodes is 
highly encouraged to constrain possible causes for the ob¬ 
served decline in activity strength. Generally speaking, 
however, 313P’s activity levels in both 2003 and 2014 
are comparable to those of other MBCs both in terms 
of peak to tal dust mas ses and peak dust-to-nucleus mass 
ratios (cf. |Hsieh||2014|). 

Our nucleus size estimate also permits us to predict 
the comet’s brightness in the absence of activity (or at 
least at the level of activity exhibited in 2004) in the 
future (Figure]^). Observed magnitudes fainter than 
predicted will indicate that the nucleus smaller than es¬ 
timated here (implying that residual activity was present 
in 2004), while magnitudes brighter than predicted will 
indicate the presence of activity. 

3.2. Morphological Analysis 

Highly variable morphology was observed for 313P in 
both 2003 and 2014. The projected orientation of the 
dust tail in the sky rotates counterclockwise during both 
active periods, each time lagging slightly behind the sim¬ 
ilarly rotating antisolar vector (as projected on the sky), 
exhibiting strong similarities t o the morphologica l evo¬ 
lution of P/La Sagra in 2010 (|Hsieh et al. 2012c). For 
P/La Sagra, two dust tails (one oriented in the antiso- 
lar direction, and another aligned with the object’s or¬ 
bit pl ane) were observed late in its observed active pe¬ 
riod (jHsieh et al.||2012c|). This distinctive mo rphology 
was also observed for 288P ( [Hsieh et al.||2012b[ ) and was 
interpreted to indicate long-lived dust emission, since 
the antisolar dust tail indicated the presence of small, 
fast-moving dust particles that must have been recently 
emitted, while the orbit-plane-aligned dust tail corre¬ 
sponded to large, slow-moving particles that must have 
been ejected much longer ago. Only one tail has been 
observed for 313P to date, but if the observed activity 
is in fact sublimation-driven, further observations could 
eventually reveal the development of a second dust tail. 

To quant itatively analyze 313P ’ s dust morphology, we 
perform a Finson & Probstein (1968) numerical dust 
modeling analysis, in this formalism, dust grains are 
released with zero velocity from the nucleus some time 
in the past, and their position is computed at the time of 
observations, accounting for solar gravity and solar radia¬ 
tion pressure. A range of ejection times and particle sizes 
are considered, where particle sizes are parametrized us¬ 
ing the ratio of solar radiation pressure to solar gravity, 
/d, where 1/P gives approximate particle radii, a, in mi¬ 
crons. The locus of the present positions of all particles 
with a common ejection time is called a synchrone, while 
the locus of the present positions of all particles with a 
common /3 is a syndyne. 

On each image, no gap is observed between the dust 
tail and the nucleus, implying either that (1) a contin¬ 
uum of dust grain sizes up to very large grains is present, 
where in practice, the seeing disk prevents us from set¬ 
ting an upper limit to the grain size, or (2) there has 
been continued emission of dust from some time in the 


past until the present, where the seeing disk does not 
allow us determine whether the emission has very re¬ 
cently stopped. The envelope of the coma is effectively 
enclosed by a syndyne representing /dmax (the minimum 
particle size) and a synchrone representing tej,min (the 
earliest ejection time), such that all the visible dust has 
P < /^max and was ejected at tej > tej,min- The values 
of /^max and tej,min themsclvcs represent lower and upper 
limits, respectively, since observations with higher S/N 
could reveal fainter surface brightness features, thus ex¬ 
tending the envelope of the observed coma and tail. 

In Figure]^ we overplot syndynes (black lines) repre¬ 
senting particles sizes ranging from P = 0.005 to /3 = 0.03 
(a ^ 30 — 200 jim) and synchrones (red lines) represent¬ 
ing ejection times from tej = —10 days to tej = —130 
days (with respect to the time of observations) on im¬ 
ages from Figure 12 We find good agreement among our 
models of individual observations indicating /dmax >0.03 
(equivalent to particles tens of microns in radius, and 
larger), and emission start times of no later than 2003 
July 24T10 {u < 10° ± 3°) for the 2003 activity, and 
no later than 2014 July 28±10 {n < 351° ± 3°) for the 
2014 activity, i.e., with activity persisting for >3 months 
each time. We find no evidence of significant changes 
(i.e., outbursts or quiescent periods) in the intensity of 
the activity or in the dust particle size distribution over 
the emission period. 

3.3. Dynamical Analysis 

To ascertain some of 313P’s basic dynami¬ 
cal properties, we first numerically propagate 
its orbit for 10 Myr using the OrbitO integrator 
(http://hamilton.dm.unipi.it/astdys/), accounting 
for the four outer planets (Jupiter to Neptune). We 
calculate proper orbital elements of Op = 3.15204 AU, 
Cp = 0.20560, and ip = 10.30°, and a Lyapunov time of 
Ti = 12 000 yr. This short Ti indicates that 313P’s orbit 
is intrinsically chaotic, perhaps due to the proximity 
of two 3-body mean motion resonances (MMRs) with 
Jupiter and Saturn, 11J-1S-5A and 10J+12S-7A. These 
two resonances are relatively weak, but have overlapping 
zones of influence for Cp > 0.15, forming a highly chaotic 
zone around 313P. 

Next, we extend our integrations to 50 Myr, aiming to 
better understand 313P’s long-term dynamical behavior. 
We find that the object’s orbital elements remain nearly 
constant over even this longer integration period (Fig¬ 
ure]^). Thus, despite residing in a chaotic region of the 
asteroid belt, the orbit of 313P appears relatively stable. 
Such behavior has als o been observed for ot her asteroids 
in high-order MMRs (Tsiganis et al. 2 002L illustrating 


a phe nomenon known as stable chao^lVlilani ^ Nobili 

We also search for an associated dynamical family to 
check whether 313P may have originate d in a recent frag¬ 
mentation event like 133P and 288P (Nesvorny et al. 
2008 Novakovic et al. |2012|). Employing the hierarchi¬ 
cal clustering method (^ppala et al.|1990), we find that 
313P is associated with t he ^155 Myr old Lixia ohna as¬ 
teroid family (Figure [Novakovic et al.pOlOj ) , joining 
the family at the low cut-off velocity of d ^ 20 m s“^, 
ma king it the second MBC in this family after P/2012 


TI (Hsieh et al. 2013). All Lixiaohna members, inclnd- 


ing 313P, have similarly chaotic orbits (Novakovic et al 
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Fig. 3. — Syndynes (black lines) and synchrones (red lines) overplotted on data from (a) 2003 September 30, (b) 2003 October 
23, (c) 2003 November 28, (d) 2014 October 1, (e) 2014 October 30, and (f) 2014 November 18, where /3 = 0.03 syndynes and 
tej = —10 day synchrones are the upper-most or left-most black and red lines, respectively, in each panel. 



Time [Myr] 



3.11 3.12 3.13 3.14 3.15 3.16 3.17 3.18 

Proper semi-major axis, ap (au] 


Fig. 4.— (a) Evolution of mean semimajor axis (top) and mean eccentricity (bottom) of 313P over 50 Myr from numerical 
integrations, (b) Plots of proper semimajor axis versus proper eccentricity (top) and sine of proper inclination (bottom) for 
Lixiaohua family members (black dots) and MBCs 313P (blue X symbol) and P/2012 T1 (green X symbol). Vertical red shaded 
regions mark the 11J-1S-5A and 10J+12S-7A MMRs. 


2010|), and as such, 313P’s Lixiaohua family association 
should be considered as reliable as those of other mem¬ 
bers, despite the chaoticity of its orbit inferred here. Un¬ 
fortunately, as in the case of P/2012 Tl, any search for 
associations with younger families will be complicated 
by the high number density of asteroids in this region 
of orbital element space, as well as the region’s intrinsic 
chaoticity. 


4. DISCUSSION 


After 133P (IHsieh et al.|2004l 

2010 

iJewitt et al.|2014b) 

and 238P ( 

Hsieh et al.| 2Ullb 

), 313 

P is now the third 


MBC observed to exhibit recurrent activity. We con¬ 
clude from this evidence, along with photometric and 
dust modeling results presented here, that 313P’s activ¬ 
ity is likely sublimation-driven, perhaps due to the col- 


lisional excavation of subsurface ice (Hsieh et al. 2004 
C apria et al.||2012] e.g .,), consistent with the conclusions 
otijewitt et 

Preliminary results from numerical dynamical model¬ 
ing indicate that some main-belt objects with moderately 
large eccentricities and inclinations similar to those of 
313P, particularly those near mean-motion resonances, 
may be temporarily-capt ured JFC-like interlopers, even 
if they have Tj > 3 (|Hsieh & Haghighipourj |2014|). 
Therefore, while our integrations of 313P itself indicate 
that it is moderately stable, and thus possibly native to 
the main belt, we caution that its dynamical origin can¬ 
not be considered absolutely certain. 

Continued observations of 313P are encouraged until 
March 2015 (^e end of the 2014-2015 observing win¬ 
dow; Figure to continue monitoring its photomet- 
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ric and morphological evolution. From September 2015 
to May 2016, 313P will have a true anomaly range of 
93° < v < 135°, and may provide opportunities to char¬ 
acterize its inactive nucleus. However, given aforemen¬ 
tioned examples of MBCs rema inin g active over simi¬ 
lar true anomaly ranges (Section 3.1), observations from 
November 2016 to July 2017 (154'' < v < 184°), may 
prove more reliable for nucleus characterization purposes. 

The heretofore unrecognized existence of clearly 
cometary images of 313P in SDSS data, despite at least 
one d edicated effort to sear ch for active comets in that 
data (Solontoi et al. 2010|), suggests that renewed ef¬ 
forts to search tor IVlBCs (and comets in general) in 
SDSS data, taking advantage of recent reprocessing of 
that data including improved blended source discrimina¬ 
tion (Solontoi 2014, private communication), could prove 
fruitful. Searches for MBCs in other data archives, espe¬ 
cially those obtained by large telescopes or with signifi¬ 
cant sky coverage, could be similarly productive. 
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